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ABSTRACT 

Spectroscopy at 8-13pm with T-ReCS on Gemini-S is presented for 3 galaxies 
with substantial silicate absorption features, NGC 3094, NGC 7172 and NGC 5506. 
In the galaxies with the deepest absorption bands, the silicate profile towards the 
nuclei is well represented by the emissivity function derived from the circumstellar 
emission from the red supergiant, \i Cephei which is also representative of the mid- 
infrared absorption in the diffuse interstellar medium in the Galaxy. There is spectral 
structure near 11.2 pm in NGC 3094 which may be due to a component of crystalline 
silicates. In NGC 5506, the depth of the silicate absorption increases from north to 
south across the nucleus, suggestive of a dusty structure on scales of 10s of parsecs. 
We discuss the profile of the silicate absorption band towards galaxy nuclei and the 
relationship between the 9.7 /im silicate and 3.4 fim hydrocarbon absorption bands. 
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1 INTRODUCTION 

Observations of active galaxies have shown that there are 
, very large columns of material along the line of sight to 
many nuclei. The standard picture is that Seyferts of type 

2 suffer significantly greater obscuration than those of type 
1, much of which probably arises in dusty circumnuclear 
structures with differing inclinations to the line of sight (eg. 
Antonucci 1993). These columns can be investigated via the 
cut-offs at X-ray wavelengths, the extinction in the visible 
and near-infrared, and dust absorption bands in the near 
and mid-infrared, all of which probe to different depths. X- 
ray observations have shown that at least 50% of nearby 
Seyfert 2 galaxies are obscured by column densities greater 
than 10 24 cm" 2 (Risaliti, Maiolini, & Salvati 1999) and that 
this material is concentrated into relatively small volumes 
(< a few parsecs). Mid-infrared observations have demon- 
strated that substantial silicate absorption depths at 9.7 pm 
are relatively common towards Seyfert 2 galaxies, but rare in 
Seyfert of type 1 (Aitken & Roche 1985, Roche et al 1991), 
again suggesting large obscuring columns towards the un- 
derlying warm emitting dust. 

Dust around AGN is likely to be affected by the hard 
photon flux from the active nucleus; for example Aitken & 
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Roche (1985) have argued that the absence of PAH emis- 
sion bands in active nuclei can be explained by destruction 
of the small grains that carry the bands. On the other hand, 
the shape of the silicate absorption feature towards the most 
heavily extinguished galaxy known, NGC 4418 (Roche et al 
1986), is similar to that produced by the diffuse interstellar 
medium in our Milky Way galaxy (Roche & Aitken 1984, 
1985, Chiar & Tielens 2006). A weak absorption band at 
3.4pm, similar to that seen in the Galactic diffuse ISM, has 
been detected towards a number of nuclei (Imanishi 2000, 
2003, Mason et al 2004), again suggesting that the dust to- 
wards those objects is similar to that in the Galaxy. Ob- 
servations with Spitzer have recently detected absorption 
features from crystalline silicates towards several deeply- 
embedded galaxies and additional absorption bands from 
hydrocarbons at 6.85 and 7.25 pm and a number of molec- 
ular species (Spoon et al 2006). 



Here we present mid-infrared spectroscopy of NGC 
3094, NGC 7172 and NGC 5506, three galaxies with signifi- 
cant silicate absorption features, and investigate the nature 
of the silicate absorption towards active nuclei at subarcsec- 
ond resolution. 
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Figure 1. Acquisition images of NGC 3094, NGC 5506 and NGC 7172 at 12 fim. The contours are linear, but the flux scale is arbitrary 
and positions are relative to the peak of the mid-IR flux. Sky conditions were unstable during observations of NGC3094, and hence the 
small scale structure of this source my not be dependable. 



2 OBSERVATIONS 

Long slit spectra between 8 and 13 |xm were obtained at 
the 8-m Gemini South telescope in clearing skies with the 
facility mid-infrared imager/spectrometer, T-ReCS (Telesco 
et al 1998), in May 2004 under programme GS-2004A-C- 
2. The seeing was variable at the time of observations, so 
the resolution achieved in the acquisition images is not well 
defined, particularly for NGC 3094 which was observed as 
the skies cleared and conditions were changing rapidly. Non- 
theless, we show these short exposures in Figure 1. We do 
not have accurate astrometric reference positions for the T- 
ReCS data and so adopt the bright core as the central ref- 
erence position. The bright nuclei were centred in the slit 
before taking the spectroscopic observations. 

Spectra were obtained with standard chop and nod 
techniques (chop throw 15 arcsee) after centering the com- 
pact nucleus in the 0.36 arcsee wide slit. The instrument was 
configured with the low resolution (11 line/mm) grating giv- 
ing a dispersion of 0.0223 (Am/pixel and a spectral resolution 
of 0.08(i.m. T-ReCS has a detector scale of 0.089 arcsee per 
pixel which provides a 25 arcsee long slit in the spatial direc- 
tion and coverage of the full N photometric band, limited by 
the N filter bandpass, in the dispersion direction. Series of A 
BAB nod positions were accumulated, with total on-source 
integrations of 10, 10 and 5 minutes on source for NGC 3094, 
NGC 5506 and NGC 7172 respectively. The chop direction 
was perpendicular to the slit, and so only the signal from 
the guided beam position was accumulated. 

Correction for telluric atmospheric absorption is with 
respect to ot Cen whilst wavelength calibration is with re- 
spect to the planetary nebula NGC 6302, which has a 
rich emission line spectrum (see Casassus, Roche & Barlow 
(2000) for the spectrum and accurate wavelengths). 

The spectra were straightened, registered in the disper- 
sion direction with reference to the sharp structure in the 
absorption band due to atmospheric ozone at 9.6 microns, 
and the dispersion calibration from NGC 6302 was applied. 
The spectra were corrected for atmospheric absorption by 
straightening and growing the alpha Cen spectrum along 
the slit, dividing it into the galaxy spectra and multiply- 
ing by a blackbody at a temperature of 5770 K (Cohen et 



al 1996). We extracted and smoothed the residual flux in 
regions ±5-10 arcsee along the slit from the spectra of the 
nuclei; this is at a very low level and presumably arises from 
beam imbalances between the chopped and nodded beams. 
The deep silicate absorption bands in NGC 3094 and NGC 
7172 give rise to very low flux levels in the 9 to 10 \im region, 
and the shape of the minimum can be distorted by this resid- 
ual background. We therefore subtracted the residual signal 
from the galaxy spectra. This only has a significant effect on 
the flux between 9.2 and 10u.m, but does appear to improve 
the match to the expected silicate profiles. 

The cores of these galaxies are fairly compact and the 
slit losses for NGC 5506 and NGC 7172 are probably similar 
to those for the reference star. Nonetheless, the calibration 
uncertainty is likely to be significant and we estimate ~20%. 

In order to search for spectral variations along the slit, 
spectra were extracted in a number of different spatial bins. 
For NGC 3094 and NGC 7172, the extracted spectra off 
the nuclei were of low signal to noise, and no significant 
differences from the on-nuclei spectra were detected. NGC 
5506 is significantly brighter, and its spectra do show spatial 
structure. Spectra were extracted by summing detector rows 
together, and are shown in figure 3 at the position of the 
central peak, at ±0.3 arcsee and at ±0.6 arcsee The noise 
in the spectra increases between 9 and 10 u.m because of 
low atmospheric transmission due to absorption by ozone in 
the upper atmosphere, while the fluxes at these wavelengths 
are decreased by silicate absorption in the source. These 
effects become increasingly apparent in the spatial positions 
furthest from the nucleus. 



3 TARGETS 

The galaxies observed all show significant silicate absorp- 
tion bands in published mid-infrared spectra with no evi- 
dence of the narrow emission bands attributed to emission 
by large polycyclic aromatic hydrocarbon molecules (PAHs) 
that trace nuclear star-forming regions (Roche et al 1991), 
and show other evidence of substantial obscuration. They 
are described briefly here. 
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Figure 2. The 8-13 \xm spectra of NGC 3094 (left) and NGC 7172 (right). The solid lines represent fits to the spectra with the Trapezium 



silicate curve while the dashed lines use that obtained from the circumstcllar shell of yi Ccp. Flux is in units of 10 



W cm 



The error bars plotted are estimated from the scatter of the points in narrow wavelength bins, the spectrum of NGC 7172 has been 
smoothed to improve the S/N ratio per point. 



3.1 NGC 3094 



3.2 NGC 5506 



NGC 3094 is a nearby (z~0.008) SBa spiral galaxy hosting a 
Scyfert 2 nucleus, and inclined by ~49° to the line of sight. 
Near-infrared images show a bright bar with a compact nu- 
cleus extended over ~1.7 arcsec (Zenner & Lenzen 1993) 
while the radio emission is fairly diffuse over 15 arcsec (Con- 
don et al 1990). The T-ReCS acquisition image indicates 
that the mid-infrared source may be extended by about 1 
arcsec at a position angle of ~ 170°, but the variable seeing 
conditions make this rather uncertain. However, the spatial 
profile within the slit, which was at position angle 0° for the 
spectroscopic observations also suggests a source extended 
on a scale of ~1 arcsec Despite the narrow slit, the T-ReCS 
spectrum has flux levels within 25% of those obtained within 
a 4 arcsec aperture by Roche et al., who detected a deep sili- 
cate absorption band with rg.7 Hm ~4.3. The 3 /im spectrum 
of NGC 3094 has been presented by Imanishi (2000) who 
found a weak 3.4 u.m hydrocarbon absorption band with a 
peak optical depth, T3.4^m ~ 0.04. Imanishi also claimed to 
detect a weak 3.3 urn PAH emission band; we see no evi- 
dence for the corresponding PAH emission features in the 
8-13 u.m spectrum. 



NGC 5506 is a Sa spiral galaxy hosting an intermediate 
Seyfert 1.9 nucleus at a redshift of 0.0062. The galaxy has 
a prominent dust lane and is almost edge-on, with an incli- 
nation angle approaching 90°. The nucleus is very compact 
in the mid-infrared, and has an optical depth T9.7^ m ~ 1.4 
and a prominent [SIV] emission line at 10.5 (im (Roche et 
al 1991, Siebenmorgen et al 2004). The nucleus is classed 
as a Seyfert 1.9 as broad wings are detected on the infrared 
lines (Blanco, Ward & Wright 1990, Nagar et al 2002). The 
X-ray spectrum suggests substantial obscuration, with a col- 
umn of material corresponding to ~ 3.4 x 10 22 cm -2 (Risaliti 
et al 2000). The hydocarbon absorption band at 3.4 u.m is 
weak with an optical depth T3.4^ m ~ 0.04 (Imanishi 2000). 
The TReCS acquisition image shows a bright compact core 
but the spectrum obtained at PA 0° has a flux level about 
50% lower than detected by Roche et al (1991) in a 4 arc- 
sec aperture, consistent with a compact core within a more 
extended emitting region. The [SIV] line emission is more ex- 
tended than the dust continuum emission along PA 0° , cor- 
responding approximately to the ionization cone axis found 
by Wilson, Baldwin & Ulvestad (1985). 
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Figure 3. The spectra of NGC 5506. Left: the spectrum of the peak (top) and at positions 0.3 and 0.6 arcsec south. Right the spectrum 
at the peak and at positions 0.3 and 0.6 arcsec north of the centre. Flux is in units of 10~ 20 Wcm~ 2 um~ 1 and the dashed and solid 
lines represent fits as described for figure 2. The spectra at 0.6 arcsec S and N have been smoothed to improve the S/N ratio per point, 
except in the vicinity of the [S IV] emission line. 



3.3 NGC 7172 

NGC 7172 is a type 2 Seyfert in a nearly edge-on Sa spiral 
galaxy at a redshift of 0.0087 (Sharpies et al 1984). It dis- 
plays a prominent dust lane running almost east-west and 
the Spitzer IRAC image shows mid-infrared emission ex- 
tending over several arcseconds in the same direction; there 
is a suggestion of extension along this same axis in the T- 
ReCS acquisition image. The nucleus of NGC 7172 is heavily 
obscured in the X-ray region with N(H) ~ 8.6 x 10 22 cm~ 2 , 
and in the mid-infrared with a silicate absorption band with 
79.7Hm ~ 2.8 (Roche et al 1991). The hydrocarbon absorp- 
tion band at 3.4 it has an optical depth T3.4 Hm ~ 0.09 (Iman- 
ishi 2000). The T-ReCS image indicates that the core is ex- 
tended on a scale of ~1 arcsec at a position angle of 90°. The 
flux captured by the T-ReCS slit at a PA of 60° is about 
one third that detected by Roche et al in their 4 arscec 
aperture, consistent with a compact core within the more 
extended emitting source detected by Spitzer. 



3.4 Fits to the spectra 

In order to quantify the behaviour of the mid-IR absorption, 
we have fit the spectra with different emission and absorp- 
tion components. The fits cannot be physically realistic in 



this complex region, but do allow us to extract some quan- 
titative data. 

Following Aitken & Roche (1982) we have fit the spec- 
tra with emission spectra suffering extinction by cool silicate 
grains. In these galaxies, inclusion of a PAH emission compo- 
nent does not improve the fits as indicated by the x 2 values 
and so the results presented employ emission and absorp- 
tion spectra consisting of black-body and silicate grains. We 
used silicate grain profiles derived from the Trapezium in 
Orion and the M supergiant p Cephei, which are taken to 
be representative of molecular clouds and the diffuse Galac- 
tic interstellar medium respectively (Roche & Aitken 1984, 
Whittet et al 1988). The goodness of fit is estimated from re- 
duced \ 2 values, with the errors estimated from the scatter 
in the data points. These allow us to differentiate between 
different combinations, but the range of values giving ade- 
quate fits can be quite large for some spectra. In particular, 
the depth of the silicate absorption depends critically on the 
assumed spectral properties of the underlying emission. Here 
we assume that the underlying emission can be represented 
by grains with a blackbody emission spectrum. 

For the most heavily obscured objects, NGC 3094 and 
NGC 7172, the reduced x values indicate that a n Cep-like 
emissivity provides a much better fit than a Trapezium-like 
emissivity profile. The results are listed in Table 2 and in- 
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Table 1. Absorption features in nearby galaxies 
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spection of Fig 2 indicates that the /i Cep curve provides a 
qualitatively much better fit to the curvature in the 8-9/im 
region than the Trapezium curve. This is in line with previ- 
ous conclusions for deeply-embedded objects, e.g. for NGC 
4418 (Roche et al 1986) or IRAS 08572+3915 (Dudley & 
Wynn- Williams 1997) and from Spitzer results (Spoon et al 
2006). The dust towards the deeply embedded galaxies is 
well-represented by dust very similar to that in the diffuse 
ISM in the Galaxy (Roche & Aitken 1985, Chiar & Tie- 
lens 2006). However, the fit to NGC 3094 shows additional 
structure near 11/xm which may represent a contribution 
from crystalline grains. Spoon et al (2006) have detected 
absorption features at 11, 16, 19, 23 and 28/im which they 
attribute to crystalline silicates in the Spitzer IRS spectra of 
galaxies with deep silicate absorption bands with T9.7 Hm > 
2.9. Unfortunately, the S/N ratio of the present data are 
insufficient to determine whether the additional 11/im crys- 
talline absorption is distributed across the emitting region 
or whether it is centrally concentrated. 

For NGC 5506, the silicate absorption feature is much 
shallower with T9.7 Hm ~1.4 with the Trapezium profile or 
1.1 with the fi Cep profile towards the nucleus. As found for 
several other galaxies with moderate silicate absorption (e.g. 
in the Circinus galaxy; Roche et al 2006), the Trapezium 
curve provides the better fit. 



4 DISCUSSION 

4.1 Spectral variations in NGC 5506 

The dust continuum emission in NGC 5506 is compact but 
sufficiently bright to allow us to measure spectral differences 
in the line emission and silicate absorption depth over the 
central 1.2 arcsec. At the distance implied by the redshift 
of NGC 5506, this corresponds to a total distance of about 
170pc. The spectra to the north have shallower silicate ab- 
sorption than on the nucleus with 79.7 Hm ~ 1.25 at 0.6 acsec 
North, while the optical depth is higher at T9.7 Hm ~ 1.7 
at 0.6 arcsec south using the Trapezium profile; the colour 
temperature of the underlying emission is ~300 K. The [SIV] 
emission line intensity decreases by a factor 3 at 0.6 arcsec 
north and by a factor 10 at 0.6 arcsec south compared to 
the central position, reflecting the changes in silicate absorp- 
tion depth and the corresponding extinction. The continuum 



omission decreases more quickly than the line emission with 
distance from the nucleus so that the equivalent width of 
the [SIV] line increases. The [SIV] intensity is lower on the 
south side of the nucleus, consistent with greater extinction 
than to the north if the intrinsic distribution is symmetric. 
The total [SIV] flux detected in the TReCS slit, 1.7 10~ 20 
W cm -2 , is 30 percent of that detected by ISO in a large 
(~ 14 x 20 arcsec) aperture (Sturm et al 2002) and about 
half that detected by Siebenmorgen et al (2004) in a 3 arc- 
sec slit. The 12.8 [im [Ne II] line is not detected in any of 
the TReCS spectra, suggesting that the [Nell] emission line 
detected by ISO arises from more spatially extended, lower 
ionization gas than the [SIV] emission. 

These spatial trends are consistent with an inclined ob- 
scuring layer that increases the extinction to both dust and 
line emitting regions to the south of the nucleus and results 
in lower extinction to the north, and which extends over at 
least 80 pc. The variations in the depth of the silicate ab- 
sorption band and the [SIV] emission line are qualitatively 
similar to those found in NGC 1068 by Mason et al (2006) 
and Rhee & Larkin (2006) and in Circinus by Roche et al 
(2006) who also found variations in the silicate absorption 
depth on subarcsec scales and line emitting regions more ex- 
tended than the continuum emission. Adopting an A(V) to 
T"9.7nm ratio of 18, as found in the diffuse ISM in the Galaxy 
(Roche & Aitken 1984, Chiar & Tielens 2006) indicates that 
the corresponding visual extinctions at 0.6 acsec S, on the 
nucleus and 0.6 arcsec N are Av ~ 30 mag, 25 mag and 22 
mag respectively. 



4.2 Absorbing Columns 

The fits to the galaxy spectra hold only for simple models 
of a warm underlying emitting nucleus suffering extinction 
by cold layers. Nonetheless, the match between the silicate 
absorption profile between 8 — 9|J.m and 11.5 — 13(tm in NGC 
3094 and that of /i Cep, representing the profile in the diffuse 
ISM in the Galaxy, is remarkable. It is noteworthy that the 
galaxies with very deep silicate absorption bands have pro- 
files similar to the diffuse ISM while those with shallower sil- 
icate absorption are better matched by the molecular cloud 
(Trapezium) curve. This is illustrated in Table 1, where the 
9.7 |0.m absorption depths are listed together with the \ 2 
values obtained with the different silicate profiles for the 
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galaxies discussed here, along with other examples of nearby 
AGN from the literature. It seems that obscuring columns 
of > 3 are better fit by the fi Cep curve, suggesting 

that molecular cloud material does not produce most of the 
absorption. The origin of the difference between the Trapez- 
ium and n Cep profiles is not clear but it may result from 
mantles on grains or, more likely, grain size effects (e.g. see 
Bowey et al 1998). We expect AGN to preferentially destroy 
small grains (Aitken & Roche 1985, Voit 1992), and so we 
might expect the grains with Trapezium-like emissivities to 
be found closest to the AGN, and small grain destruction 
to be most complete in objects with the lowest absorbing 
columns. Alternatively, it is possible that the silicate grains 
closest to the nuclei are located in very dense clumps, where 
they are cold and can accrete molecular material and may 
coagulate, and perhaps these regions are seen preferentially 
in the less obscured nuclei. 

NGC 5506 and NGC 7172 are close to edge-on to the 
line of sight, with prominent dust lanes in optical images, so 
it is perhaps likely that some of this extinction arises in the 
galactic ISM. The columns of material implied from the X- 
ray absorption (~ 3.4 x 10 22 and ~ 8.6 x 10 22 cm~ 2 ; Risaliti 
ct al 2000) are similar to those inferred from the silicate 
absorption bands towards the nuclei (~ 2.8 x 10 22 and ~ 
6.4 x 10 22 cm~ 2 for NGC 5506 and NGC 7172 respectively, 
estimated from the relationship between A(V) and rg.7^ m 
in the Galactic ISM (Roche & Aitken 1984) and a ratio 
of Nh / Ay = 1.9 10 21 cm" 2 mag" 1 (Bohlin, Savage & 
Drake 1978). The similar absorbing columns suggest that the 
mid-IR and X-ray emitting regions suffer similar extinctions 
in these objects, unlike many other AGN where the X-ray 
column is much higher than estimates from the IR. Roche & 
Aitken (1985) have found that the ratio of silicate absorption 
depth to visual extinction is higher towards the Galactic 
centre than towards more local objects in the Milky way, so 
the comparisons in absorbing column derived from different 
wavebands are uncertain by at least a factor 2. 

The 3.4/im absorption band is attributed to aliphatic 
hydrocarbons and is detected along sightlines through the 
diffuse interstellar medium rather than through molecular 
material (e.g. Chiar et al 2000). The ratio of TgAp.m to T9.7^ m 
provides estimates of the relative columns of hydrocarbons 
to silicates. Along the line of sight to the Galactic centre and 
more local sightlines within the Galaxy this ratio is ~ 0.06 
and does not show significant variations (Chiar et al 2000). 
The ratio in the galaxies listed in table 1 varies between 0.01 
and 0.17. If the 3.4/um band is representative of diffuse ISM 
material, we might expect the ratio of T3.4 Hm to r^.r^m to 
be closer to the Galactic value in the galaxies that show the 
narrow /i Cep like silicate profile, but this is not the case. In 
the most heavily obscured galaxies, the 3.4 (im hydrocarbon 
absorption band is weaker relative to the silicate absorp- 
tion band than in our Galaxy (see also Imanishi 2000). This 
may reflect differences in grain populations, perhaps as a 
result of irradiation by the AGN. Spectropolarimetry shows 
that the grains producing the 3.4 fim absorption towards 
the Galactic centre are unpolarized in contrast to those pro- 
ducing the silicate absorption (Chiar et al 2006, Aitken et 
al 1986), indicating that the hydrocarbon and silicate ab- 
sorption bands arise from separate populations. However, 
contributions from different components, differences in the 



emitting volumes at 3 and 10 /tm and complex circumnu- 
clear structures may also be important. 



5 CONCLUSIONS 

NGC 5506 shows variations in silicate absorption depth on a 
scales of tens of parsecs. This could arise in a dusty structure 
extending over many tens of parsecs in the nuclear region. 
High spatial resolution observations of NGC 1068 and Circi- 
nus show similar effects, so that such structures appear to 
be a common property of circumnuclear material in nearby 
AGN. 

The 10 /im silicate absorption in galaxies with deep 
absorption closely resembles the silicate profile in the dif- 
fuse interstellar medium in our Galaxy, as represented by 
the circumstellar emission from the red supergiant fi Cep 
while less heavily obscured nuclei appear to be better fit by 
the Trapezium emissivity curve, which is representative of 
silicates in molecular clouds. This could reflect preferential 
destruction of small grains in the galaxies with the lowest 
absorbing columns. However, the absorption depth of the 
3.4/im aliphatic hydrocarbon band, which is prominent in 
the Galactic diffuse ISM is not well correlated with the sili- 
cate absorption depth, even in galaxies with diffuse-ISM-like 
profiles. 

In common with many heavily obscured ULIRGS, NGC 
3094 shows evidence of spectral structure at 11 /im which 
may arise from a population of crystalline silicates; higher 
quality spectra may be able to determine whether these 
grains are confined to the nuclear region or are more widely 
distributed. 
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